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1  INTRODUCTION 


Over  the  past  twenty  years,  it  has  become  clear  that  hmits  on  per¬ 
formance  of  low-frequency  sonar-systems  from  oceanographic  variabihty  are 
measurable,  and  that  internal  waves  are  a  major  contributor  to  these  limits. 
However,  in  deep  water  these  hmits  have  been  so  weak  that  sonar  systems 
in  the  field  have  not  been  built  of  a  size  and  type  that  would  be  affected  by 
decoherence  due  to  internal  waves. 

Two  changes  in  the  situation  suggest  a  re-investigation  of  the  limits 
due  to  internal  waves.  First,  target  submarines  have  gotten  quieter,  so  that 
larger  and  more  ambitious  sonar  systems  with  higher  gain  are  being  contem¬ 
plated.  Second,  shallow-water,  coastal  regions  are  taking  on  a  higher  level 
of  importance  in  operational  considerations.  Both  of  these  changes  move 
internal-wave  effects  into  a  more  significant  position. 

This  report  will  present  the  state  of  knowledge  of  internal-wave  effects 
from  the  point  of  view  of  low-frequency  sonar  systems,  and  will  make  rec¬ 
ommendations  to  the  sonar-system  designer  and  operator  as  to  when  to  ex¬ 
pect  internal-wave  effects  to  become  important.  It  will  also  explain  where 
oceanographic  measurement  and  analysis  are  needed  to  determine  the  hmits 
on  sonar-system  design  and  performance  that  are  due  to  internal  waves. 

The  calculations  in  this  report  are  based  on  the  imiversal  Garrett-Munk 
spectrum,  with  differences  in  different  parts  of  the  ocean  arising  from  the  be¬ 
havior  of  different  water  masses.  Although  different  strengths  of  the  internal- 
wave  field  are  easily  taken  into  account,  if  the  internal- wave  field  is  signifi¬ 
cantly  different  from  the  GM  model,  these  calculations  will  not  be  accurate. 
This  particularly  may  be  the  case  in  some  shallow-water  locations.  Further¬ 
more,  other  effects  such  as  solitons,  solibores,  or  intrusions  may  be  important 
in  some  locations.  Nevertheless,  we  feel  that  it  is  worthwhile  doing  the  GM 
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calculation  as  a  reference  from  which  to  judge  any  alternative  models  that 
may  arise  in  the  future. 

The  effectiveness  of  internal  waves  in  causing  acoustic  fluctuations  is  de¬ 
pendent  on  oceanographic  behavior  (the  temperature-salinity  (T-S)  relation) 
whose  systematics  would  need  to  be  established  more  precisely  than  is  known 
at  present.  It  is  also  dependent  on  the  local  strength  of  internal  waves,  whose 
variability  is  known  to  exist,  but  whose  systematics  need  to  be  established 
also.  Finally  it  is  dependent  on  the  local  sound-speed  profile,  which  needs  to 
be  locally  measured,  and  the  buoyancy-frequency  profile,  which  needs  to  be 
at  least  modelled  reasonably  precisely. 

A  Scenario 

Imagine  that  a  sonar  system  is  built  that  has  10^  hydrophones  in  it,  and 
that  it  is  deployed  near  a  coast  in  1000  meters  of  water.  Without  oceano¬ 
graphic  variability,  such  a  system  is  capable  of  50  dB  of  signal  gain  against 
isotropic  noise.  At  100  Hz  it  would  have  an  aperture  of  1000  m  in  the  ver¬ 
tical  and  about  6  km  in  the  horizontal.  The  calculations  in  this  report  will 
show  that  at  a  range  of  500  km,  a  reference  internal-wave  model  will  pre¬ 
dict  acoustic  coherence  lengths  that  vary  widely  between  ocean  regions,  but 
with  typical  values  of  5000  m  in  the  vertical  and  25  km  in  the  horizontal. 
Thus  it  would  appear  that  a  sonar-system  designer  is  justified  in  ignoring 
internal-wave  effects. 

However,  in  our  scenario,  we  happen  to  be  in  the  North  Polar  Pacific 
(near  Japan),  and  a  nearby  storm  has  been  raging  for  a  week  or  so.  As  a  re¬ 
sult,  we  are  in  an  acoustically  sensitive  region,  and  the  internal-wave  strength 
has  gone  up  by  a  factor  of  ten  in  energy  (displacements  have  increased  by 
a  factor  of  three).  The  internal- wave-induced  coherence  lengths  drop  by  a 
factor  of  ten  to  500  m  and  2.5  km,  and  the  gain  drops  below  the  ideal  50  dB 
by  about  6  dB  (at  500  km).  In  addition,  the  sonar  system  has  been  relying 
on  some  incoherent  time  integration,  and  the  gain  available  has  dropped  also. 
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As  a  result,  a  quiet  target  has  dropped  below  threshold  and  is  now  invisible. 


As  a  modification  of  our  scenario,  we  move  in  toward  shallower  water, 
say  200-m  depth,  but  we  consider  ranges  of  20  km  instead  of  500  km.  Internal 
waves  will  not  limit  the  vertical  coherence,  because  the  water  depth  is  so 
small;  we  assume  coherence  up  to  the  maximum  of  200  m  in  the  vertical. 
Horizontal  coherence  lengths  in  the  standard  model  are  in  the  hundreds  of 
kilometers,  so  that  no  internal-wave  effects  are  significant. 

This  kind  of  analysis  is  the  starting  point  for  evaluating  the  actual  limits 
imposed  by  internal  waves.  It  seems  clear  that  internal-wave  effects  should 
be  considered  first  in  the  long-range,  deep-water  scenario,  and  that  shallow- 
water  scenarios  at  low  frequency  that  have  internal-wave  effects  will  be  much 
rarer,  involving  special  cases.  The  long-range  signal-gain  limits  due  to  in¬ 
ternal  waves  are  going  to  be  in  the  regime  of  50  dB,  with  a  more  exact 
prediction  requiring  the  monitoring  of  internal-wave  strength  as  a  function 
of  geographical  location  and  time,  in  addition  to  already-realized  monitor¬ 
ing  requirements  on  the  sound-speed  and  buoyancy-frequency  profiles.  Thus, 
knowledge  of  the  effectiveness  of  a  sonar  system  that  is  large  enough  to  be 
affected  by  internal  waves  will  require: 


•  monitoring  the  local  sound-speed  profile  (already  done). 

•  knowing  the  local  T-S  relation  in  an  acoustically  relevant  manner. 

•  monitoring  the  local  buoyancy-frequency  profile. 

•  monitoring  the  local  strength  of  internal  waves. 

The  rest  of  this  report  will  describe  our  present  knowledge  of  all  these 
factors,  and  how  they  enter  into  sonar-system  performance. 


2  OCEAN  CHARACTERIZATION 


2.1  Sound  Speed  and  Buoyancy  Frequency 


In  order  to  make  a  realistic  prediction  of  acoustic  behavior,  including  the 
effects  of  internal  waves,  it  is  necessary  to  use  realistic  profiles  with  depth  of 
the  sound  speed  and  buoyancy  frequency.  Such  profiles  can  be  obtained  from 
a  well-known  database  that  has  been  compiled  under  the  name  of  Levitus. 
This  database  has  information  on  temperature  and  salinity  as  a  function  of 
depth  for  each  season  and  on  a  one-degree  grid  of  latitude  and  longitude. 
Information  is  provided  every  ten  meters  in  depth  down  to  1500  m.  There  is 
also  a  yearly  average  database,  called  “annual” ,  which  includes  information 
down  to  the  depth  of  the  ocean,  which  is  in  most  places  between  4000  and 
6000  m. 

In  order  to  develop  realistic  profiles  that  represent  a  large  ocean  region, 
such  as  the  North  Pacific,  we  can  pick  a  particular  latitude  (say  SON)  and 
average  the  profiles  over  longitude  (say  140W-200W).  However,  it  turns  out 
that  this  is  a  dangerous  and  unreliable  procedure,  because  the  act  of  averag¬ 
ing  the  sound-speed  profile  and  the  buoyancy  profile  separately  can  generate 
two  profiles  that  are  unrepresentative  in  important  ways.  Therefore  it  is  bet¬ 
ter  to  look  at  the  profiles  of  the  different  longitudes  separately  and  pick  a 
particular  longitude  that  is  “representative” . 

The  representative  geographical  locations  we  have  picked  for  each  of 
a  number  of  ocean  regions  are  shown  in  Table  1.  The  sound-speed  and 
buoyancy  profiles  from  the  Levitus  database  that  correspond  to  these  repre¬ 
sentative  locations  are  shown  in  Figure  1. 
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Table  1:  Ocean  regions  chosen  for  this  study. 


Ocean  Region 

Latitude 

Longitude 

North  Polar  Atlantic 

60°iV 

40W 

North  Polar  Pacific 

50°N 

m°w 

North  Atlantic 

30°A 

30^W 

North  Pacific 

30°N 

175° VP 

Equatorial  Atlantic 

0° 

25°W 

Equatorial  Pacific 

O'’ 

150° VP 

Equatorial  Indian 

0° 

90°E 

South  Atlantic 

30°5 

15°IP 

South  Pacific 

30°S 

135°W 

South  Indian 

30°S 

75°E 

South  Polar  Atlantic 

60°S 

20°W^ 

South  Polar  Indian 

60°S 

120°E 

The  sound-speed  profile  will  control  the  deterministic  part  of  the  behav¬ 
ior  of  acoustic  propagation  in  the  ocean  region  being  considered.  The  North 
Atlantic  has  a  considerably  weaker  thermocline  than  the  other  regions,  and 
of  course  the  polar  regions  have  almost  no  thermocline,  instead  maintaining 
a  reverse  gradient  for  depths  almost  all  the  way  to  the  surface. 

The  buoyancy-frequency  profile  will  be  needed  for  predicting  the  effects 
of  internal  waves.  It  varies  in  the  upper  ocean  by  factors  of  order  five  between 
different  ocean  basins,  and  this  factor  will  have  dramatic  effects  on  internal- 
wave  acoustic  effects. 


2.2  Internal- Wave  Spectrum  and  Strength  ((^) 


Deep-water  internal  waves  are  characterized  reasonably  well  by  their 
power  spectrum,  which  has  been  empirically  determined  by  Garrett  and 
Munk  (GM)  to  be  expressible  as  factorizable  in  frequency  and  vertical  mode 
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Depth 


Ocean— Basin  Levitus  annual  data 


Sound— speed  (m/s)  Buoyancy  Frequency(cph) 


Figure  1:  The  profiles  of  sound  speed  and  buoyancy  frequency,  for  the  ocean 
regions  used  in  this  study,  from  the  Levitus  annual  database.  The  various 
profiles  are  listed  from  longest  to  shortest  dash  (solid  line  first);  North  At¬ 
lantic,  South  Atlantic,  North  Pacific,  South  Pacific,  South  Indian,  North 
Polar  Atlantic,  South  Polar  Atlantic,  South  Polar  Indian. 
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number:  [1] 


S{j,co)  =  ieiz)) 


P  +  P 


(2-1) 


where  Nj  and  are  normalization  constants  defined  in  Flatte  and  Esswein,  [2] 
such  that  the  sum  over  j  and  integral  over  u  of  the  spectrum  gives 


Later,  the  dispersion  relation  of  internal  waves  relating  wavevector  k  to 
mode  number  j  and  frequency  lo  will  be  needed.  Define  k  as  the  magnitude 
of  the  horizontal  component  of  the  wavevector:  k  =  ^k^  -|-  k^.  The  vertical 
wavevector  component  is  k^  and  the  frequency  is  co. 


k  =  aj{u‘^ ky  =  ajn^.  (2-2) 


The  quantity  a  depends  on  the  vertical  profile  of  buoyancy  frequency  n{z), 
and  is  approximately  given  by: 


a 


-1 


1 

TT 


(2-3) 


where  D  is  the  depth  of  the  ocean.  This  integral  of  the  buoyancy-frequency 
profile  will  be  an  important  parameter  in  determining  the  acoustic  effects  of 
internal  waves.  The  value  of  or'^  for  the  deep  ocean  is  typically  between  3 
and  5  cph-km,  or  in  other  units,  5  to  9  m  s“L  Its  value  in  shallow  water  can 
be  estimated  by  using  a  depth  of  200  m  and  a  buoyancy  frequency  of  8  cph, 
which  gives  1.6  cph-km,  or  3  m  s“^.  Thus  internal  waves  in  200  m  of  water 
move  at  speeds  that  typically  are  only  2  or  3  times  slower  than  the  speed  of 
internal  waves  in  deep  water. 


The  quantity  in  Equation  (2-1)  is  the  displacement  of  internal  waves  at 
a  particular  depth;  it  is  related  to  the  energy  density  of  internal  waves  and  to 
the  buoyancy  profile.  It  has  been  determined  by  measurements  that  (C^('2)) 
is  near  50  m^  at  depths  in  the  deep  ocean  at  which  the  buoyancy  frequency 
is  3  cph.  Let  Co  be  this  reference  level,  with  an  accompanying  definition  of 
no  as  3  cph.  However,  there  are  substantial  variations  from  this  “reference” 
level.  These  variations  are  of  great  interest  to  oceanographers,  because  they 


8 


represent  the  sources  and  sinks  of  the  internal  wave  field.  They  are  also 
of  great  interest  to  the  sonar-system  designer,  because  the  acoustic  effects 
of  internal  waves  will  be  proportional  (in  some  sense)  to  the  internal-wave 
strength. 


2.2.1  The  Equator  Problem 


Here  we  mention  a  difficulty  in  the  definition  of  the  GM  spectrum. 
The  expressions  for  most  quantities  depend  on  the  inertial  frequency,  which 
goes  to  zero  at  the  equator.  Therefore,  the  frequency  dependence  of  the 
internal- wave  displacement  spectrum  is  singular  at  the  equator.  There  is 
no  agreed-upon  simple  modification  of  the  GM  spectrum  that  avoids  this 
problem,  although  the  physics  of  the  various  kinds  of  waves  that  are  involved 
(Kelvin,  Rossby,  etc)  is  well  known.  Until  a  simple  GM  modification  for 
equatorial  regions  is  available,  the  sonar-system  designer  has  little  guidance 
in  this  area. 

2.3  Ratio  of  Sound-Speed  to  Buoyancy  Effects  (F) 


Even  if  the  internal- wave  reference  displacement  Co  applied  universally 
to  the  world’s  oceans,  the  acoustic  effects  due  to  internal  waves  would  vary 
from  ocean  basin  to  ocean  basin.  This  is  because  different  ocean  basins 
have  different  water  masses.  That  is,  different  ocean  basins  have  different 
temperature-salinity  relations.  Acoustic  effects  of  internal  waves  depend  on 
the  product  of  potential  sound-speed  gradient  and  internal-wave  displace¬ 
ment,  while  internal- wave  strength  depends  on  the  product  of  density  gradi¬ 
ent  and  internal- wave  displacement.  Let  us  call  F  the  ratio  of  acoustic  effect 
to  internal-wave  effect: 

^  AC,/C, 

ApIp 
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where  Cp  is  the  potential  sound  speed  (that  is,  the  sound  speed  with  the 
pressure  effect  removed).  Both  potential  sound  speed  and  density  depend 
on  temperature,  salinity,  and  pressure,  according  to  the  equation  of  state 
of  seawater.  There  are  well-known  formulas  for  these  quantities,  based  on 
extensive  measurements  in  the  laboratory:  notably,  the  results  of  Chen  and 
Millero  [3]  for  density  and  Del  Grosso  [4]  for  sound  speed.  In  order  to  use 
these  results,  the  temperature-saUnity  characteristics  of  ocean  basins  must 
be  oceanographically  understood. 


2.3.1  Water  Masses 


Different  ocean  regions  have  different  temperature-salinity  (TS)  rela¬ 
tions,  which  to  first  order  are  independent  of  time  (including  season).  Figure 
2  shows  the  TS  relations  for  twelve  ocean  basins  by  plotting  S  vs  T  for  a 
particular  latitude  and  ocean  basin,  with  each  plot  having  several  different 
longitudes.  It  is  seen  that  the  different  longitudes  fall  closely  on  each  other 
over  a  substantial  depth  region.  Although  it  isn’t  shown  in  the  plots,  the 
depth  below  which  the  different  longitudes  are  quite  similar  is  about  200  me¬ 
ters.  Above  that  depth,  the  seasonal  heating  of  the  sun  imposes  substantial 
differences. 

If  the  plots  of  salinity  or  temperature  were  presented  as  functions  of 
depth,  there  would  be  much  more  difference  between  different  longitudes. 
By  presenting  the  results  in  terms  of  S  and  T,  the  oceanographers  have 
thus  been  able  to  identify  (by  ocean  basin  and  depth  region)  different  water 
masses  that  change  very  little  over  the  years.  Each  of  the  curves  of  salinity 
versus  temperature  can  be  used  to  determine  the  F  factor  as  a  function  of 
temperature. 
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A  summary  plot  of  TS  relations  for  different  ocean  basins  is  shown 
in  Figure  3,  where  representative  longitudes  from  the  previous  figures  have 
been  chosen. 


2.3.2  Results  for  F 


It  is  tempting  to  start  by  approximating  the  equation  of  state  of  seawa¬ 
ter  by  first  derivative  terms,  [5]  but  that  turns  out  to  be  inadequate  in  this 
situation.  (The  points  in  the  Levitus  database  are  too  far  apart  to  use  first 
derivatives,  and  the  coefficients  of  the  first-derivative  terms  are  too  variable 
with  temperature  and  depth.)  For  the  results  shown  here,  the  complete  for¬ 
mulas  have  been  used.  Figure  4  shows  F  as  a  function  of  depth  for  a  number 
of  ocean  regions.  It  is  seen  that  F  varies  over  a  factor  of  three,  and  that 
will  have  dramatic  effects  on  the  relative  acoustic  effects  of  internal  waves 
in  different  oceans.  Because  different  seasons  and  latitudes  have  different 
temperatures  at  the  same  depth,  the  value  of  F  varies  over  longitude  and 
season.  However,  if  F  is  plotted  as  a  function  of  temperature,  the  underlying 
water-mass  characteristic  of  the  ocean  basin  is  revealed  more  clearly  by  a 
narrowing  of  the  F  curves.  Figure  5  shows  F  as  a  function  of  temperature 
for  a  number  of  ocean  regions,  and  Figure  6  shows  F  for  representative  lon¬ 
gitudes  from  each  of  those  ocean  regions.  Again  we  see  the  variation  of  F 
over  a  factor  of  three. 


2.4  Shallow- Water  Modifications 


Internal  waves  in  shallow  water  are  the  subject  of  on-going  research. 
However,  for  the  sonar-system  designer,  and  the  oceanographic-experiment 
designer,  it  is  important  to  have  a  qualitative  picture  of  the  shallow-water 
situation.  Measmrements  of  internal  waves  in  water  with  depths  down  to 
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Salinity 


Figure  4:  Acoustic-eflFect  factor  F  as  a  function  of  depth  for  twelve  ocean 
regions.  This  factor  is  the  ratio  of  fractional  sound-speed  variation  to  the 
fractional  density  variation  due  to  an  internal  wave.  The  different  curves  on 
each  plot  are  for  different  longitudes. 
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Figure  5:  Acoustic-effect  fax:tor  F  as  a  function  of  temperatm:e  for  twelve 
ocean  regions.  This  factor  is  the  ratio  of  fractional  sound-speed  variation  to 
the  fractional  density  variation  due  to  an  internal  wave.  The  different  curves 
on  each  plot  are  for  different  longitudes. 
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Ocean— Basin  Levitus  annual  data 
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Figure  6:  Acoustic-effect  factor  F  as  a  function  of  temperature  for  the  ocean 
regions  used  in  this  study.  The  various  profiles  are  listed  from  longest  to 
shortest  dash  (solid  line  first):  North  Atlantic,  South  Atlantic,  North  Pacific, 
South  Pacific,  South  Indian,  North  Polar  Atlantic,  South  Polar  Atlantic, 
South  Polar  Indian. 
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100  meters  have  been  made  in  a  variety  of  areas;  those  measurements  have 
shown  that  the  GM  spectrum  is  not  an  unreasonable  model  when  a  few 
modifications  are  taken  into  account.  (The  most  recent  measmements  have 
been  made  off  the  East  and  West  coasts  of  the  United  States.) 


•  In  water  of  a  few  hundred  meters  depth  the  internal  waves  are  vertically 
coherent  to  a  high  degree.  This  means  that  the  vertical  correlation 
length  of  internal  waves  is  effectively  the  water  depth.  But  in  deep 
water  the  internal-wave  vertical  correlation  length  is  about  200  meters. 
Since  this  is  not  far  from  the  depth  in  shallow  water,  the  GM  model  as 
it  stands  is  not  inappropriate. 

•  Observations  in  shallow  water  show  internal-wave  rms  displacements 
of  order  10  meters,  just  as  in  deep  water.  However,  the  buoyancy  fre¬ 
quency  in  shallow  water  is  systematically  larger  than  in  the  surface 
waters  of  the  deep  ocean  by  about  a  faotor  of  two  (the  difference  be¬ 
tween,  say,  8  and  4  cph).  This  does  not  change  the  spectrum,  but  does 
affect  the  overall  strength  of  internal-wave  acoustic  effects  in  a  simple 
way. 

•  Shallow-water  observations  have  identified  internal  waves  other  than 
the  superposition  of  linear  waves  described  by  the  GM  spectrum.  In 
particular,  solitons  and  solibores  are  often  observed.  [6,  7]  This  study 
will  not  evaluate  the  effect  of  these  other  internal  waves.  However, 
the  SWARM  experiment  shows  an  example  in  which  the  GM  part  of 
the  spectrum  has  about  the  same  rms  displacement  as  the  solitons, 
and  since  solitons  have  a  much  shorter  correlation  length  along  their 
direction  of  propagation,  their  effects  on  acoustics  are  often  much  less. 
When  they  are  weak,  as  seems  to  be  the  case  on  the  Pacific  coast 
(Murray  Levine,  private  communication,  1997),  they  have  negligible 
effect.  When  they  are  strong,  (e.g.  in  the  SWARM  data)  their  effects 
would  be  no  more  than  comparable  to  GM  effects  in  size,  though  their 
frequencies  are  much  higher  than  GM  internal  waves.  The  sonar-system 
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designer  is  justified  for  the  most  part  in  keeping  them  in  abeyance  for 
the  moment.  When  future  research  establishes  systematic  models  for 
solitons  and  solibores  that  are  comparable  to  the  GM  model,  then 
sonar-system  designers  will  be  able  to  incorporate  their  effects. 


A  set  of  shallow-water  locations  around  the  world  where  a  combination 
of  political  and  military  factors  indicate  the  possible  existence  of  future  naval 
conflict  are  given  in  Table  2,  along  with  their  typical  ocean-depth  regimes. 

Note  that  only  the  Falklands  are  in  the  Atlantic,  or  part  of  any  “Atlantic 
Rim  !  Only  the  Falklands  and  the  related  Strait  of  Magellan  are  in  the 
Southern  Hemisphere. 


Table  2:  Interesting  shallow-water  ocean  regions. 


Ocean  Region 

Latitude 

Longitude 

Depth  Regime 

Persian  Gulf 

27°N 

32°E 

30-90  m 

Strait  of  Hormuz 

20"iV 

30-90  m 

Taiwan  Strait 

25°iV 

120"E 

50-80  m 

Strait  of  Malacca 

5°N 

100"£^ 

40-1000  m 

Molucca  Sea 

0° 

125"F; 

1000-4000  m 

Korea  Strait 

35°iV 

130".E 

50-200  m 

Panama  Canal 

10°N 

80W 

30-1000  m 

Strait  of  Magellan 

53°S 

14°W 

30-500  m 

The  Falklands 

52°S 

Q2°W 

60-500  m 

2.5  Sound-Speed- Variance  Depth  Dependence 


The  effects  of  sound-speed  variability  in  ocean  “surface”  waters  (i.e. 
those  waters  within  about  300  m  of  the  surface)  are  important.  The  GM 
model  was  put  together  to  apply  to  the  “deep”  ocean:  i.e.  waters  at  depths 
of  1000  m  or  more.  At  these  large  depths  it  is  very  useful  to  apply  the  WKBJ 
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approximation  to  internal-wave  vertical  wavefunctions.  Yet  the  GM  model 
has  been  remarkably  useful  for  depths  much  closer  to  the  surface.  If  it  is 
applied  all  the  way  to  the  surface,  the  WKBJ  approximate  wavefunctions 
do  not  go  to  zero  as  the  real  internal- wave  wavefunction  must.  Therefore  a 
correction  must  be  made  to  reduce  the  calculated  effects  of  internal  waves 
when  the  WKBJ  approximations  are  used.  [8] 

Figure  7  shows  the  variance  of  sound-speed  fluctuations  in  an  internal- 
wave  model  in  which  the  buoyancy  frequency  has  a  Munk  canonical  profile.  [9] 
The  WKBJ  approximation  is  seen  to  increase  from  large  depths  all  the  way 
to  the  surface.  A  numerical  simulation  of  internal  waves  with  proper  vertical 
mode  functions  results  in  a  sound-speed  variance  that  goes  to  zero  at  the 
surface,  and  an  empirical  function  m{z)  has  been  developed  to  correct  the 
WKBJ  approximation  to  account  for  the  numerical-simulation  results:  [8] 

m(z)  =  1  —  (2-5) 

where  zi  =  55  m.  This  empirical  function  has  been  shown  to  result  in  accurate 
predictions  of  internal-wave-induced  sound-speed  acoustic  fluctuations  for 
acoustic  energy  that  travels  near  the  ocean  surface  during  part  of  its  path.  [8] 
This  empirical  function  has  been  used  in  all  the  calculations  in  Section  3  of 
this  report. 

2.6  Internal- Wave  Correlation  Lengths 


In  the  ray  approximation,  acoustic  effects  depend  on  the  integral  of 
sound-speed  variations  along  the  direction  of  propagation  of  a  ray.  Upon 
evaluation  of  such  acoustic  effects  by  a  path-integral  technique  (see  Section 
3),  a  property  of  importance  appears  in  the  form  of  a  correlation  length  along 
the  ray,  which  has  been  called  Lp.[5]  For  almost  twenty  years,  formulas  have 
been  used  for  Lp  that  were  presented  by  Munk  and  Zachariasen,  based  on 
approximating  rays  locally  as  straight  lines.  [10]  Recently  a  generalization  to 
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Figure  7:  The  internal- wave-induced  sound-speed  variance  (/i^)  as  a  function 
of  depth,  for  a  case  in  which  the  profile  of  buoyancy  frequency  is  canonical. 
The  sohd  fine  is  the  analytic  profile  expected  from  an  exponential  buoy¬ 
ancy  frequency  and  the  WKBJ  approximation.  The  dashed  line  is  the  result 
of  numerical  simulation  of  an  internal-wave  field  with  use  of  the  GM  spec¬ 
trum,  and  with  accurate  mode  wave  functions  calculated  on  the  basis  of  the 
buoyancy-frequency  profile.  The  dash-dot  line  is  the  WKBJ  approximation 
corrected  with  an  empirical  function. 
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this  approximation  that  accounts  for  ray  curvature  has  been  developed.  [8] 

The  correlation  length  Lp  can  be  defined  in  terms  of  two  dimensionless 
parameters.  The  first  parameter  Ps  is  related  to  the  anisotropy  of  the  internal- 
wave  spectrum: 

Ps  =  —  tan  6.  (2-6) 

<^i 

where  Us  is  the  buoyancy  firequency  at  the  depth  at  which  Lp  is  being  eval¬ 
uated,  and  ^  is  an  angle  with  the  horizontal  at  that  depth. 


An  expression  for  Lp  is  then: 


where 


F(/?)  =  2/3  r 

Jo 


J(j{v)v^dv 


(2-7) 


(2-8) 


[v^  +  pf 

and  where  rig  is  an  “effective”  buoyancy  frequency  that  applies  to  two  points 
separated  along  an  acoustic  ray.  The  details  of  the  derivation  of  this  formula 
are  presented  in  the  reference;  here  it  suffices  to  present  the  results  as  they 
are  used  in  our  calculations,  and  to  say  that  if  Ue/ris  is  set  to  imity,  the  old 
straight-fine  approximation  is  recovered. 


Note  that  F{P)  is  defined  such  that 

/*oo 

/  F(P)d0  =  1  (2-9) 

J  0 


The  behavior  of  Ue/ris  is  controlled  by  the  depth  behavior  of  two  points 
on  a  ray  trajectory  symmetrically  placed  on  either  side  of  the  evaluation 
point,  and  by  the  depth  dependence  of  the  buoyancy  frequency.  A  series 
of  approximations,  assuming  the  ray  trajectory  to  have  constant  curvature, 
results  in  the  definition  of  another  dimensionless  parameter  cr: 
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where  B  is  a  scale  length  of  the  buoyancy-frequency  profile:  typically  1300 
m.  The  parameter  a  represents  the  ratio  between  an  internal-wave  length 
scale  (rj)  and  the  ray  radius  of  curvature  (r).  A  typical  value  of  Tj  is  1500 
km. 


The  result  of  all  these  definitions  is  a  final  expression  for  Lp  in  terms  of 
dimensionless  parameters: 


Lp 


/4Brjy^  N^Nj 
V  3  /  jj{P  +  j^) 


cos{/3ps)F{/3), 


(2-11) 


Figure  8  shows  Lp  os  a.  function  of  Ps  (corresponding  to  ray  angle  with 
the  horizontal)  evaluated  for  various  values  of  cr  (corresponding  to  curvature 
of  the  ray  trajectory).  The  a  values  axe  0,  1,  10,  20,  50,  100,  200,  and  500. 
The  values  of  Lp  for  <  1  are  reasonably  close  to  those  from  the  straight- 
line  approximation.  The  values  of  Lp  for  more  curved  rays  (larger  a)  are 
more  isotropic:  smaller  in  the  horizontal  direction  and  larger  for  steep  rays. 


The  results  shown  in  Figure  8  were  calculated  numerically  by  tabulating 
F(/3).  It  is  useful  to  find  an  analytic  expression  for  Lp  in  order  to  make 
computer  calculations  more  efficient.  A  useful  expression  has  been  found: [8] 


Lp  =  L 


l-exp[- 


(2-12) 


where  the  values  of  the  various  parameters  are:  Pc=3.5,  (7c=0.0204,  p=0.385, 
Pa=0.5,  g=1.3,  and  5a=2.0.  The  value  of  LpQ  is  given  by: 


Lpo  — 


N^N,M. 


-1 


(2-13) 


where  Nj,  Mj,  and  are  normalization  constants  defined  in  Flatte  and 
Esswein.[2] 


The  above  empirical  formula  applies  for  a  >  6.  If  a  <  6  then  the  formula 
is  used  with  a  =  6.  Figure  9  shows  results  from  the  analytic  empirical  formula 
for  cr  values  of  0,  6,  10,  20,  50,  100,  200,  and  500. 
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Values  of  a  for  ocean  rays  depend  upon  the  sound-speed  profile  and  the 
parameters  characterizing  the  internal- wave  field  (r/).  Figme  10  shows  two 
sets  of  profiles  (for  sound  speed  and  buoyancy  frequency)  for  conditions  in 
the  north  Pacific.  The  first  set  corresponds  to  a  canonical  profile  used  in 
simulation  of  the  SLICE89  experiment.  [11]  The  second  set  is  an  average  of 
profiles  taken  along  a  3.2  Mm  path  between  California  and  Hawaii  for  the 
ATOC  experiment.  [12]  Figure  11  shows  <7  as  a  function  of  depth  for  the 
two  profile  sets  in  the  north  Pacific.  The  values  of  a  range  up  to  600,  which 
corresponds  to  a  ray  radius  of  curvature  of  about  2  km. 

2.7  Geographical  and  Seasonal  Variations  -  The  Sys¬ 
tematic  View  (F  vs 


In  order  to  predict  the  effects  of  internal  waves  on  acoustic  propagation, 
a  reasonable  pair  of  sound-speed  and  buoyancy  profiles  need  to  be  obtained, 
which  we  assume  to  happen  either  from  using  a  database  like  that  of  Levitus, 
or  by  measinrement  in  the  region  and  time  of  interest.  In  addition  to  these 
two  profiles,  two  parameters  must  be  provided:  the  strength  of  the  internal- 
wave  spectrum  and  the  F  factor.  It  therefore  becomes  of  interest  to  have 
a  systematic  view  of  these  two  parameters  as  a  function  of  geography  and 
geophysical  time  (season,  month,  or  whatever  is  relevant  to  possible  changes 
in  these  two  parameters). 

We  have  very  little  information  about  variations  in  Co  21  function  of 
time  or  space.  The  best  we  can  do  at  present  is  use  one  number  (the  reference 
level)  or  refer  to  very  specific  experiments.  For  example,  the  SLICE89  exper¬ 
iment  in  the  North  Pacific  between  California  and  Hawaii  in  1989  measured 
a  Co  of  about  25  m^[13],  while  preliminary  results  from  the  ATOC  experi¬ 
ment  in  a  nearby  region  in  1997  give  about  100  m^  (Kevin  Heaney,  private 
communication,  1997).  The  GM75  reference  level  of  50  m^  was  observed  in 
a  nearby  region  in  1974  by  several  experiments.  Thus  variations  of  a  factor 
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Sound  — speed  (m/s)  Buoyancy  Frequency(cph) 


Figure  10:  The  profiles  of  sound  speed  and  buoyancy  frequency,  for  the  aver¬ 
age  of  a  set  of  profiles  observed  in  the  ocean  between  California  and  Hawaii 
(dashed),  and  for  a  canonical  profile  used  in  simulations  of  the  SLICE89 
experiment  (solid). 
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Figure  11:  The  profiles  of  a  for  the  two  sets  of  profiles  shown  in  the  previous 
figure. 
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of  four  with  different  geophysical  time  have  been  observed  in  an  ocean-basin 
sized  region.  These  experiments  varied  in  their  spatial  averaging  from  a  few 
kilometers  to  5  megameters. 

The  variation  of  F  with  time  may  be  assumed  to  be  small  for  depths 
deeper  than  300  m,  and  thus  a  single  map  of  F  can  be  created.  Figure  12 
shows  such  a  map  created  from  the  previous  plots  in  this  study.  A  temper¬ 
ature  of  10®  C  has  been  chosen  as  representative  of  depths  of  interest.  The 
value  of  F  varies  by  a  factor  of  three,  although  in  temperate  latitudes  the 
variation  is  considerably  less. 


Ocean-Basin  F  value  at  10deg  C  (except  polar  regions) 


300*  0*  60*  120*  180*  240* 


300*  0*  60*  120*  180*  240* 


Figure  12:  Geographical  distribution  of  F  at  a  temperature  of  10®  C. 
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3  ACOUSTIC  SIGNAL-GAIN  ESTIMATES 


Calculations  of  acoustic  fluctuations  have  been  carried  out  by  means  of 
formulas  derived  by  use  of  path-integral  techniques.  [5, 14, 15]  These  formulas 
comprise  numerical  integrals  along  rays  from  a  source.  These  rays  could  be 
characterized  in  various  ways;  in  this  study  we  have  used  a  source  at  500-m 
depth  in  order  to  obtain  rays  whose  upper  turning  points  (UTP)  span  the 
region  between  the  surface  and  500-m  depth.  It  has  been  found  that  the 
variation  with  UTP  depth  can  be  substantial;  however,  we  will  pick  a  single 
number  that  characterizes  an  acoustic  fluctuation  quantity  without  worrying 
too  much  about  UTP-depth  dependence  of  the  result. 

Calculations  have  been  made  for  ranges  of  500  km  and  2000  km.  The 
results  scale  with  range  according  to  relatively  simple  expectations,  so  that 
extrapolations  to  other  ranges  are  not  difficult.  Therefore  we  can  present  the 
results  only  for  500  km  and  refer  to  scaling  to  obtain  results  at  other  ranges. 

Calculations  have  been  made  for  twelve  ocean  regions.  We  have  studied 
the  equatorial  regions  as  a  function  of  latitude,  and  have  used  values  between 
5°  and  10°  to  avoid  the  singularity  at  0°. 

Five  acoustic  fluctuation  quantities  have  been  calculated.  Each  is  ex¬ 
pressed  in  a  form  that  is  independent  of  acoustic  frequency  to  first  approx¬ 
imation.  However,  the  sonar-system  designer  and  operator  are  of  course 
interested  in  dependences  on  acoustic  frequency,  which  will  be  discussed. 

Each  quantity  is  calculated  separately  for  each  ray.  The  actual  acoustic 
field  at  the  receiver  will  be  a  sum  of  rays:  deterministic  multipath.  The 
combination  of  rays  implied  by  a  given  source-receiver  geometry  is  assumed 
to  be  a  solvable  signal-processing  problem. 
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Each  quantity  and  its  implications  on  sonar-system  design  and  operation 
are  described  below: 

•  critical  frequency  fc  -  The  implication  of  r  is  mainly  in  controlling  the 
upper  limit  on  frequencies  for  which  internal-wave  effects  are  negligible. 
Taking  a  simple  rule-of-thumb  that  wavefront  fluctuations  less  than 
A/ 10  are  negligible,  the  critical  frequency  can  be  calculated  from  the 
rms  travel-time  fluctuation  r.  Thus,  fc  =  100/r  where  fc  is  in  Hz 
and  r  is  in  ms.  For  example,  if  r  =  4  ms,  then  frequencies  below  25  Hz 
have  negligible  internal-wave  effects.  Conversely,  internal-wave  effects 
for  frequencies  above  25  Hz  will  be  non-negligible. 

•  coherent  bandwidth  cbw  -  Incoherence  between  frequencies  can  be  trans¬ 
lated  into  the  time  domain  as  the  spread  of  an  individual  pulse  dmring 
transmission.  Considering  that  the  width  of  a  pulse  is  about  twice  the 
rms  spread  tq,  the  coherent  bandwidth  can  be  calculated  from  the  rms 
spread  as  cbw  =  500/ro  where  cbw  is  in  Hz  and  tq  is  in  ms.  For 
example,  a  4-ms  rms  spread  translates  into  a  coherent  bandwidth  of 
125  Hz. 

•  coherence  time  to,  expressed  in  acoustic  periods  -  The  coherence  time  is 
controlled  by  the  motion  of  internal  waves,  and  it  controls  the  maximum 
integration  time  that  can  be  used  for  coherent  temporal  integration. 
For  example,  a  coherence  time  of  40,000  periods  at  100  Hz  translates 
to  400  seconds  coherent  integration  time. 

•  vertical  coherence  length  l^,  expressed  in  acoustic  wavelengths  -  The 
vertical  coherence  length  controls  the  maximum  length  that  can  be  used 
for  coherent  spatial  integration.  In  the  vertical,  plane-wave  beamform- 
ers  are  not  adequate;  complete  models  of  the  propagation  to  form  a 
“matched  field”  must  be  done,  but  the  principle  is  the  same;  the  signal 
gain  available  from  a  vertical  aperture  is  approximately  the  vertical 
aperture  in  wavelengths,  sot  he  signal  gain  in  the  vertical  is  l^.  Note 
that  in  the  case  of  shallow  water,  the  water  depth  may  be  less  than  the 
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physical  length  represented  by  in  that  case  the  water  depth  in  units 
of  acoustic  wavelengths  replaces  ly. 

•  horizontal  coherence  length  4,  expressed  in  acoustic  wavelengths  -  The 
horizontal  coherence  length  is  in  principle  even  simpler  than  the  verti¬ 
cal,  because  plane  waves  are  a  good  approximation  in  the  horizontal. 
Thus  the  signal  gain  in  the  horizontal  is  Ih- 

A  two-dimensional,  billboard  array  that  takes  advantage  of  both  vertical 
and  horizontal  aperture  will  be  able  to  achieve  a  signal  gain  (SG)  against 
isotropic  noise  of  lylh-  The  final  result  of  our  calculations  will  be  a  plot  of 
SG  versus  range  for  a  variety  of  ocean  regions,  assuming  specific  values  of  F 
and  The  sonar-system  designer  can  use  this  information  to  incorporate 
limits  resulting  from  internal-wave  effects  into  a  design.  The  sonar-system 
operator,  with  local  information  about  F  and  Co)  can  predict  whether  an 
operating  system  will  be  impeded  strongly  by  internal  waves  or  not. 

In  particular,  given  the  same  sound-speed  profile,  each  of  the  acoustic 
fiuctuation  quantities  listed  above  are  inversely  proportional  to  a  factor  Q, 
where: 

Q  =  F-  [n{z)f  ■  Crms  (3-1) 

Note  that  the  factor  is  Crms,  not  Co-  The  local  internal- wave  displacement 
Crms  does  depend  on  n{z),  but  it  is  more  direct  to  use  the  local  displacement 
in  order  to  more  easily  use  experimental  results.  A  systematic  example  of 
this  is  the  assumption  we  will  make  that  Crms  is  about  the  same  in  shallow 
water  as  in  deep,  while  the  buoyancy  frequency  n{z)  is  typically  a  factor  of 
two  larger  in  shallow  water  than  in  deep. 

In  deep  water,  signal  gain  is  inversely  proportional  to  because  each 
of  ly  and  Ih  yield  one  power  of  Q.  In  shallow  water,  the  vertical  factor  is 
missing,  giving  only  one  power  of  Q. 
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3.1  Definition  of  Ray  and  Ray- Tube  Geometries 


It  is  useful  to  describe  the  sound  speed  in  the  ocean  as:  [5] 

c(s)  =  Co[l  +  U{z)  +  (3-2) 

where  Co  is  the  axial  (minimum)  sound  speed,  U {z)  is  a  function  of  vertical 
position  only  and  represents  the  deterministic  sound  channel,  and  ^Ji{^  is  the 
zero-mean  fractional  sound-speed  variation  caused  by  internal  waves. 

dxx\Zray{x^  —  7^^z[^(^)]  (3”3) 

The  trajectory  of  a  ray  in  the  ocean  sound  charmel,  where  2:  is  the  depth 
and  X  is  the  range,  satisfies  the  Eikonal  equation.  Rather  than  define  many 
variables,  we  give  the  parabolic  approximation  to  that  equation  with  a  range- 
independent  sound  speed  c(z),  since  it  wiU  lead  to  simpler  ra5d;ube  functions 
later: 

Because  we  are  using  a  geometrical-optics  basis  for  our  calculations,  the 
Green’s  function  takes  the  form  of  ray-tube  behavior.  [5]  Figure  13  shows 
a  curved  ray  from  a  source  to  a  receiver  with  definitions  of  two  ray-tube 
functions  and  ^2-  These  functions  of  range  between  somce  and  receiver 
begin  with  value  zero  and  end  with  value  unity.  The  difference  between  them 
is  that  begins  at  the  source  and  ^2  begins  at  the  receiver.  Each  of  the  ^ 
functions  satisfy  the  raytube  equation: 

dxx^{x)^-U"[Zray{x)\^{x)  =  0  (3-4) 

Later,  a  function  51(0;)  will  appear  in  some  of  the  expressions  for  acous¬ 
tic  fluctuation  quantities;  it  is  the  solution  of  Equation  (3-4)  with  initial 
conditions  <S'i(0)  =  0  and  <5i(a;)  =  1..  where  U"  is  the  second  derivative  of 
U(z). 
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^l(R) 

'  I 


Figure  13:  Ray  and  ra5d;ube  definitions.  The  function  starts  from  zero 
at  the  source  and  has  a  value  of  unity  at  the  receiver.  The  function  ^2  is 
zero  at  the  receiver  and  unity  at  the  source.  These  two  functions  represent 
the  difference  between  two  rays  launched  from  the  source  or  receiver  with 
infinitesimally  different  launch  angles. 


3.2  Acoustic  Fluctuation  Quantities 


In  the  following,  the  acoustic  frequency  will  be  denoted  by  cr^.  (The 
logarithmic  dependence  of  cbw,  1^,  and  Ih  on  acoustic  frequency  is  weak;  we 
have  used  100  Hz  for  a  a  to  evaluate  these  expressions.) 

Various  quantities  averaged  over  the  internal- wave  spectrum  will  be  de¬ 
noted  by  {quantity}.  Detailed  expressions  for  these  averages  are  given  in 
Flatte  and  Stoughton.  [15] 


3.2.1  /c  (critical  frequency  for  significant  internal-wave  effects) 


Expressions  for  fc  are  obtained  from  its  relation  to  r,  the  rms  travel-time 
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fluctuation: 


_  100 
Jc  — 

T 


where  fc  is  in  Hz  and  r  is  in  ms.  The  quantity  r  is  given  by: 

1  rR 

9  J-  /  ,  i  O.  ^ 


rtt 

/  dx{ix^)Lp 
J  0 


3.2.2  cbw  (coherent  bandwidth) 


The  coherent  bandwidth  is  obtained  from  its  relation  to  tq,  the  rms 
pulse  spread; 

cbw  = -  (3-7) 

where  cbw  is  in  Hz  and  tq  is  in  ms. 


The  quantity  tq  is  given  by: 

2  _  1 

°  2^  a 


I”  dx{^^)L,{kl} 
dx\{f,'f)L;{K^]Ux)Ux')f 


3.2.3  to  (coherence  time  in  periods) 


The  quantity  to  is  given  by: 


f  dx{/j,^)Lp{u‘^} 

J  0 
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3.2.4  ly  (vertical  coherence  length  in  wavelengths) 


The  quantity  ly  is  given  by: 

ly^  =  [ln(<Tar)]2  f  dx{n‘^)Lp{kl}[^i{x)f  (3-10) 

J  0 


3.2.5  Ih  (horizontal  coherence  length  in  wavelengths) 


The  quantity  4  is  given  by  a  rather  complicated  formula,  because  of  the 
complicated  nature  of  the  internal- wave  spectrum  as  a  function  of  horizontal 
wavenumber  components. [15] 

where  Ay  is  a  reference  horizontal  separation.  The  calculations  here  have 
been  done  with  Ay  set  to  1  km;  note  that  the  dependence  of  4  on  Ay  is 
quite  weak. 


3.3  Calculational  Results  for  Acoustic  Fluctuation  Quan¬ 
tities 


Table  3  shows  the  results  for  calculations  of  the  five  acoustic  fluctuation 
quantities,  as  described  above.  Some  systematic  behavior  may  be  noted: 


•  Very  different  sound-speed  profiles  in  high-latitude  areas  result  in  the 
critical  frequencies  being  significantly  higher  there. 

•  The  behavior  of  the  acoustic  propagation  in  polar  regions  where  there 
is  usually  no  sound  channel  results  in  significantly  lower  coherent  band- 
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Table  3:  Acoustic  fluctuation  quantities  for  500-km  range  through  a 
reference-level  internal-wave  field.  The  units  are:  Hz  for  the  critical  fre¬ 
quency  fc  (below  which  internal-wave  effects  are  negligible)  and  the  coherent 
bandwidth  cbw;  kiloperiods  for  the  coherence  time  to;  and  wavelengths  for 
the  coherence  lengths  ly  and  Ih- 


Ocean  Region 

— 

NBA 

100 

25 

250 

300 

14000 

NPP 

14 

10 

250 

150 

2500 

NA 

25 

500 

90 

600 

3000 

NP 

16 

250 

60 

500 

2000 

EA 

20 

500 

180 

1000 

1700 

EP 

20 

500 

200 

800 

1700 

El 

200 

600 

1700 

SA 

90 

600 

3000 

SP 

20 

170 

80 

600 

3000 

SI 

20 

700 

80 

600 

3000 

SPA 

70 

50 

100 

200 

8000 

SPI 

50 

17 

150 

200 

8000 

width  at  those  frequencies  for  which  internal-wave  effects  are  not  neg¬ 
ligible. 

•  Temperate-latitude  regions  are  reasonably  similar  to  each  other. 

•  Equatorial  regions  have  larger  coherence  times,  larger  vertical  coher¬ 
ence  lengths,  and  smaller  horizontal  coherence  lengths  than  nearby 
temperate  regions. 


3.4  Billbocird  Signal  Gains  for  Different  Conditions 
{SG  vs  R) 


By  use  of  the  values  of  Ih  and  ly  from  Table  3,  the  signal  gain  (SG) 
against  isotropic  noise  can  be  calculated.  The  coherence  lengths  to  a  good 
approximation  scale  as  the  square  root  of  the  range.  The  signal  gain  is  the 
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product  of  two  coherence  lengths,  so  it  scales  linearly  with  range.  This  scaling 
is  reasonable  as  long  as  one  remains  in  deep  water.  As  one  moves  into  shallow 
water  two  modifications  need  to  be  made. 


•  The  value  of  C  will  be  assumed  to  stay  constant,  while  the  value  of  n{z) 
increases  by  a  factor  of  two.  By  use  of  Equation  (3-1),  the  horizontal 
coherence  length  will  decrease  by  a  factor  of  four  below  its  value  ex¬ 
trapolated  from  deep  water.  This  will  decrease  the  signal  gain  by  6 
dB. 

•  The  vertical  coherence  length  extrapolated  from  deep  water  will  be 
replaced  by  the  water  depth.  To  represent  shallow  water,  a  depth  of 
200  m  will  be  used  to  calculate  SG.  This  gives  a  different  correction  for 
each  ocean  region,  but  the  typical  correction  is  a  reduction  in  SG  of  20 
dB. 


Values  of  SG  for  different  ocean  regions  as  a  function  of  range  are  shown 
in  Figure  14.  The  shallow-water  point  has  been  chosen  to  be  20  km,  and  the 
ranges  between  20  and  500  km  have  been  filled  in  by  a  smooth  curve,  not  by 
any  precise  calculation. 

Achievement  of  these  maximal  signal  gains  would  require  the  develop¬ 
ment  of  billboard  acoustic  arrays  that  are  tens  of  kilometers  wide  and  span 
the  ocean  water  column  in  depth.  Hundreds  of  thousands  of  hydrophones  in 
each  array  would  be  necessary.  Nevertheless,  these  are  the  ocean  limits. 
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Ocean— basin  annual  data 
F  =  20  f  =  7  m 
Source  and  Receiver  are  both 
R  from  shore  and  R  from  each  other 


10’  102  103 

Range  (km) 

Figure  14:  Limits  due  to  a  reference-level  internal- wave  field  (C^  =  50  m^;  F 
=  20)  for  signal  gain  against  isotropic  noise  for  a  billboard  array.  Vertical  and 
horizontal  apertures  are  set  to  the  coherence  lengths  due  to  internal  waves. 
Results  are  shown  for  different  ocean  regions  as  a  function  of  range  between 
source  and  receiver.  The  short  ranges  are  assumed  to  be  in  shallow  water  as 
explained  in  the  text.  (In  deep  water  the  maximum  signal  gain  is  inversely 
proportional  to  range.) 
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4  RECOMMENDATIONS 


This  study  has  combined  present  oceanographic  knowledge  of  the  tempera¬ 
ture-salinity  structure  in  the  sea  with  knowledge  of  the  effects  of  internal 
waves  on  acoustic  propagation.  It  has  been  revealed  that  two  parameters 
are  of  paramount  importance  to  the  understanding  of  internal-wave  effects 
on  sonar  systems:  the  local  strength  of  the  internal-wave  field  and  the  lo¬ 
cal  value  of  a  factor  F  that  represents  the  ratio  of  fractional  variations  in 
sound-speed  and  density.  The  F  factor  has  been  oceanographically  studied 
in  the  form  of  temperature  sahnity  relationships,  and  reasonable  values  for 
many,  but  not  all,  ocean  regions  are  available.  The  variation  of  internal- 
wave  strength  is  known  to  be  substantial,  but  no  systematic  information  is 
available  at  the  present  time  on  which  to  base  designs  of  sonar  systems. 

If  the  behavior  of  sonar  systems  in  conditions  in  which  they  are  limited 
by  internal  waves  is  desired,  the  following  recommendations  are  strongly 
suggested: 


•  Initiate  a  substantial  monitoring  program  to  follow  the  strength  of  in¬ 
ternal  waves  as  a  function  of  season  and  ocean  region. 

•  Validate  and  make  more  precise  the  use  of  the  F  factor  described  in 
this  report. 

•  Develop  an  appropriate  generalization  of  the  GM  spectrum  that  is  valid 
across  the  equator. 

•  Verify  the  validity  of  the  use  of  the  GM  spectrum  in  shallow  water,  or 
replace  it  with  a  new  systematic  model. 
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